The hazards of genotype imputation when mapping disease susceptibility variants.

BACKGROUND: The cost-free increase in statistical power of using imputation to infer missing genotypes is undoubtedly appealing, but is it hazard-free? This case study of three type-2 diabetes (T2D) loci demonstrates that it is not; it sheds light on why this is so and raises concerns as to the shortcomings of imputation at disease loci, where haplotypes differ between cases and reference panel. RESULTS: T2D-associated variants were previously identified using targeted sequencing. We removed these significantly associated SNPs and used neighbouring SNPs to infer them by imputation. We compared imputed with observed genotypes, examined the altered pattern of T2D-SNP association, and investigated the cause of imputation errors by studying haplotype structure. Most T2D variants were incorrectly imputed with a low density of scaffold SNPs, but the majority failed to impute even at high density, despite obtaining high certainty scores. Missing and discordant imputation errors, which were observed disproportionately for the risk alleles, produced monomorphic genotype calls or false-negative associations. We show that haplotypes carrying risk alleles are considerably more common in the T2D cases than the reference panel, for all loci. CONCLUSIONS: Imputation is not a panacea for fine mapping, nor for meta-analysing multiple GWAS based on different arrays and different populations. A total of 80% of the SNPs we have tested are not included in array platforms, explaining why these and other such associated variants may previously have been missed. Regardless of the choice of software and reference haplotypes, imputation drives genotype inference towards the reference panel, introducing errors at disease loci.

The hazards of genotype imputation in chromosomal regions under selection: A case study using the Lactase gene region.

Although imputation of missing SNP results has been widely used in genetic studies, claims about the quality and usefulness of imputation have outnumbered the few studies that have questioned its limitations. But it is becoming clear that these limitations are real-for example, disease association signals can be missed in regions of LD breakdown. Here, as a case study, using the chromosomal region of the well-known lactase gene, LCT, we address the issue of imputation in the context of variants that have become frequent in a limited number of modern population groups only recently, due to selection. We study SNPs in a 500 bp region covering the enhancer of LCT, and compare imputed genotypes with directly genotyped data. We examine the haplotype pairs of all individuals with discrepant and missing genotypes. We highlight the nonrandom nature of the allelic errors and show that most incorrect imputations and missing data result from long haplotypes that are evolutionarily closely related to those carrying the derived alleles, while some relate to rare and recombinant haplotypes. We conclude that bias of incorrectly imputed and missing genotypes can decrease the accuracy of imputed results substantially.

New Insights Into Mitochondrial Dysfunction at Disease Susceptibility Loci in the Development of Type 2 Diabetes.

This study investigated the potential genetic mechanisms which underlie adipose tissue mitochondrial dysfunction in Type 2 diabetes (T2D), by systematically identifying nuclear-encoded mitochondrial genes (NEMGs) among the genes regulated by T2D-associated genetic loci. The target genes of these 'disease loci' were identified by mapping genetic loci associated with both disease and gene expression levels (expression quantitative trait loci, eQTL) using high resolution genetic maps, with independent estimates co-locating to within a small genetic distance. These co-locating signals were defined as T2D-eQTL and the target genes as T2D cis-genes. In total, 763 cis-genes were associated with T2D-eQTL, of which 50 were NEMGs. Independent gene expression datasets for T2D and insulin resistant cases and controls confirmed that the cis-genes and cis-NEMGs were enriched for differential expression in cases, providing independent validation that genetic maps can identify informative functional genes. Two additional results were consistent with a potential role of T2D-eQTL in regulating the 50 identified cis-NEMGs in the context of T2D risk: (1) the 50 cis-NEMGs showed greater differential expression compared to other NEMGs and (2) other NEMGs showed a trend towards significantly decreased expression if their expression levels correlated more highly with the subset of 50 cis-NEMGs. These 50 cis-NEMGs, which are differentially expressed and associated with mapped T2D disease loci, encode proteins acting within key mitochondrial pathways, including some of current therapeutic interest such as the metabolism of branched-chain amino acids, GABA and biotin.

A novel human pain insensitivity disorder caused by a point mutation in ZFHX2.

Chronic pain is a major global public health issue causing a severe impact on both the quality of life for sufferers and the wider economy. Despite the significant clinical burden, little progress has been made in terms of therapeutic development. A unique approach to identifying new human-validated analgesic drug targets is to study rare families with inherited pain insensitivity. Here we have analysed an otherwise normal family where six affected individuals display a pain insensitive phenotype that is characterized by hyposensitivity to noxious heat and painless bone fractures. This autosomal dominant disorder is found in three generations and is not associated with a peripheral neuropathy. A novel point mutation in ZFHX2, encoding a putative transcription factor expressed in small diameter sensory neurons, was identified by whole exome sequencing that segregates with the pain insensitivity. The mutation is predicted to change an evolutionarily highly conserved arginine residue 1913 to a lysine within a homeodomain. Bacterial artificial chromosome (BAC) transgenic mice bearing the orthologous murine p.R1907K mutation, as well as Zfhx2 null mutant mice, have significant deficits in pain sensitivity. Gene expression analyses in dorsal root ganglia from mutant and wild-type mice show altered expression of genes implicated in peripheral pain mechanisms. The ZFHX2 variant and downstream regulated genes associated with a human pain-insensitive phenotype are therefore potential novel targets for the development of new analgesic drugs.awx326media15680039660001.

World-wide distributions of lactase persistence alleles and the complex effects of recombination and selection.

The genetic trait of lactase persistence (LP) is associated with at least five independent functional single nucleotide variants in a regulatory region about 14 kb upstream of the lactase gene [-13910*T (rs4988235), -13907*G (rs41525747), -13915*G (rs41380347), -14009*G (rs869051967) and -14010*C (rs145946881)]. These alleles have been inferred to have spread recently and present-day frequencies have been attributed to positive selection for the ability of adult humans to digest lactose without risk of symptoms of lactose intolerance. One of the inferential approaches used to estimate the level of past selection has been to determine the extent of haplotype homozygosity (EHH) of the sequence surrounding the SNP of interest. We report here new data on the frequencies of the known LP alleles in the 'Old World' and their haplotype lineages. We examine and confirm EHH of each of the LP alleles in relation to their distinct lineages, but also show marked EHH for one of the older haplotypes that does not carry any of the five LP alleles. The region of EHH of this (B) haplotype exactly coincides with a region of suppressed recombination that is detectable in families as well as in population data, and the results show how such suppression may have exaggerated haplotype-based measures of past selection.

High-Resolution Genetic Maps Identify Multiple Type 2 Diabetes Loci at Regulatory Hotspots in African Americans and Europeans.

Interpretation of results from genome-wide association studies for T2D is challenging. Only very few loci have been replicated in African ancestry populations and the identification of the implicated functional genes remain largely undefined. We used genetic maps that capture detailed linkage disequilibrium information in European and African Americans and applied these to large T2D case-control samples in order to estimate locations for putative functional variants in both populations. Replicated T2D locations were tested for evidence of being regulatory hotspots using adipose expression. We validated a sample of our co-location intervals using next generation sequencing and functional annotation, including enhancers, transcription, and chromatin modifications. We identified 111 additional disease-susceptibility locations, 93 of which are cosmopolitan and 18 of which are European specific. We show that many previously known signals are also risk loci in African Americans. The majority of the disease locations appear to confer risk of T2D via the regulation of expression levels for a large number (266) of cis-regulated genes, the majority of which are not the nearest genes to the disease loci. Sequencing three cosmopolitan locations provided candidate functional variants that precisely co-locate with cell-specific chromatin domains and pancreatic islet enhancers. These variants have large effect sizes and are common across populations. Results show that disease-associated loci in different populations, gene expression, and cell-specific regulatory annotation can be effectively integrated by localizing these effects on high-resolution genetic maps. The cis-regulated genes provide insights into the complex molecular pathways involved and can be used as targets for sequencing and functional molecular studies.

The genomic and phenotypic diversity of Schizosaccharomyces pombe.

Natural variation within species reveals aspects of genome evolution and function. The fission yeast Schizosaccharomyces pombe is an important model for eukaryotic biology, but researchers typically use one standard laboratory strain. To extend the usefulness of this model, we surveyed the genomic and phenotypic variation in 161 natural isolates. We sequenced the genomes of all strains, finding moderate genetic diversity (π = 3 × 10(-3) substitutions/site) and weak global population structure. We estimate that dispersal of S. pombe began during human antiquity (∼340 BCE), and ancestors of these strains reached the Americas at ∼1623 CE. We quantified 74 traits, finding substantial heritable phenotypic diversity. We conducted 223 genome-wide association studies, with 89 traits showing at least one association. The most significant variant for each trait explained 22% of the phenotypic variance on average, with indels having larger effects than SNPs. This analysis represents a rich resource to examine genotype-phenotype relationships in a tractable model.

ABCC5 transporter is a novel type 2 diabetes susceptibility gene in European and African American populations.

Numerous functional studies have implicated PARL in relation to type 2 diabetes (T2D). We hypothesised that conflicting human association studies may be due to neighbouring causal variants being in linkage disequilibrium (LD) with PARL. We conducted a comprehensive candidate gene study of the extended LD genomic region that includes PARL and transporter ABCC5 using three data sets (two European and one African American), in relation to healthy glycaemic variation, visceral fat accumulation and T2D disease. We observed no evidence for previously reported T2D association with Val262Leu or PARL using array and fine-map genomic and expression data. By contrast, we observed strong evidence of T2D association with ABCC5 (intron 26) for European and African American samples (P = 3E-07) and with ABCC5 adipose expression in Europeans [odds ratio (OR) = 3.8, P = 2E-04]. The genomic location estimate for the ABCC5 functional variant, associated with all phenotypes and expression data (P = 1E-11), was identical for all samples (at Chr3q 185,136 kb B36), indicating that the risk variant is an expression quantitative trait locus (eQTL) with increased expression conferring risk of disease. That the association with T2D is observed in populations of disparate ancestry suggests the variant is a ubiquitous risk factor for T2D.

Refinement in localization and identification of gene regions associated with Crohn disease.

The risk of Crohn disease (CD) has a large genetic component. A recent meta-analysis of 6 genome-wide association studies reported 71 chromosomal intervals but does not account for all of the known genetic contribution. Here, we refine localization of the previously reported intervals and also identify additional CD susceptibility genes using a mapping approach that localizes causal variants based on genetic maps in linkage disequilibrium units (LDU maps). Using 2 of the 6 cohorts, 66 of the 71 previously reported loci are confirmed and more precise location estimates for these intervals are given. We identify 78 additional gene regions that pass genome-wide significance, providing strong evidence for 144 genes. Additionally, 56 nominally significant signals, but with more stringent and precise colocalization, are identified. In total, we provide evidence for 200 gene regions confirming that CD is truly multifactorial and complex in nature. Many identified genes have functions that are compatible with involvement in immune/inflammatory processes and seem to have a large effect in individuals with extra ileal as well as ileal inflammation. The precise locations and the evidence that some genes reflect phenotypic subgroups will help identify functional variants and will lead to greater insight of CD etiology.

A candidate-gene association study for berry colour and anthocyanin content in Vitis vinifera L.

Anthocyanin content is a trait of major interest in Vitis vinifera L. These compounds affect grape and wine quality, and have beneficial effects on human health. A candidate-gene approach was used to identify genetic variants associated with anthocyanin content in grape berries. A total of 445 polymorphisms were identified in 5 genes encoding transcription factors and 10 genes involved in either the biosynthetic pathway or transport of anthocyanins. A total of 124 SNPs were selected to examine association with a wide range of phenotypes based on RP-HPLC analysis and visual characterization. The phenotypes were total skin anthocyanin (TSA) concentration but also specific types of anthocyanins and relative abundance. The visual assessment was based on OIV (Organisation Internationale de la Vigne et du Vin) descriptors for berry and skin colour. The genes encoding the transcription factors MYB11, MYBCC and MYC(B) were significantly associated with TSA concentration. UFGT and MRP were associated with several different types of anthocyanins. Skin and pulp colour were associated with nine genes (MYB11, MYBCC, MYC(B), UFGT, MRP, DFR, LDOX, CHI and GST). Pulp colour was associated with a similar group of 11 genes (MYB11, MYBCC, MYC(B), MYC(A), UFGT, MRP, GST, DFR, LDOX, CHI and CHS(A)). Statistical interactions were observed between SNPs within the transcription factors MYB11, MYBCC and MYC(B). SNPs within LDOX interacted with MYB11 and MYC(B), while SNPs within CHI interacted with MYB11 only. Together, these findings suggest the involvement of these genes in anthocyanin content and on the regulation of anthocyanin biosynthesis. This work forms a benchmark for replication and functional studies.

Comparative analysis of genome-wide association studies signals for lipids, diabetes, and coronary heart disease: Cardiovascular Biomarker Genetics Collaboration.

AIMS: To evaluate the associations of emergent genome-wide-association study-derived coronary heart disease (CHD)-associated single nucleotide polymorphisms (SNPs) with established and emerging risk factors, and the association of genome-wide-association study-derived lipid-associated SNPs with other risk factors and CHD events. METHODS AND RESULTS: Using two case-control studies, three cross-sectional, and seven prospective studies with up to 25 000 individuals and 5794 CHD events we evaluated associations of 34 genome-wide-association study-identified SNPs with CHD risk and 16 CHD-associated risk factors or biomarkers. The Ch9p21 SNPs rs1333049 (OR 1.17; 95% confidence limits 1.11-1.24) and rs10757274 (OR 1.17; 1.09-1.26), MIA3 rs17465637 (OR 1.10; 1.04-1.15), Ch2q36 rs2943634 (OR 1.08; 1.03-1.14), APC rs383830 (OR 1.10; 1.02, 1.18), MTHFD1L rs6922269 (OR 1.10; 1.03, 1.16), CXCL12 rs501120 (OR 1.12; 1.04, 1.20), and SMAD3 rs17228212 (OR 1.11; 1.05, 1.17) were all associated with CHD risk, but not with the CHD biomarkers and risk factors measured. Among the 20 blood lipid-related SNPs, LPL rs17411031 was associated with a lower risk of CHD (OR 0.91; 0.84-0.97), an increase in Apolipoprotein AI and HDL-cholesterol, and reduced triglycerides. SORT1 rs599839 was associated with CHD risk (OR 1.20; 1.15-1.26) as well as total- and LDL-cholesterol, and apolipoprotein B. ANGPTL3 rs12042319 was associated with CHD risk (OR 1.11; 1.03, 1.19), total- and LDL-cholesterol, triglycerides, and interleukin-6. CONCLUSION: Several SNPs predicting CHD events appear to involve pathways not currently indexed by the established or emerging risk factors; others involved changes in blood lipids including triglycerides or HDL-cholesterol as well as LDL-cholesterol. The overlapping association of SNPs with multiple risk factors and biomarkers supports the existence of shared points of regulation for these phenotypes.

Dissecting the genetics of complex inheritance: linkage disequilibrium mapping provides insight into Crohn disease.

Family studies for Crohn disease (CD) report extensive linkage on chromosome 16q and pinpoint NOD2 as a possible causative locus. However, linkage is also observed in families that do not bear the most frequent NOD2 causative mutations, but no other signals on 16q have been found so far in published genome-wide association studies. Our aim is to identify this missing genetic contribution. We apply a powerful genetic mapping approach to the Wellcome Trust Case-Control Consortium and the National Institute of Diabetes and Digestive and Kidney Diseases genome-wide association data on CD. This method takes into account the underlying structure of linkage disequilibrium (LD) by using genetic distances from LD maps and provides a location for the causal agent. We find genetic heterogeneity within the NOD2 locus and also show an independent and unsuspected involvement of the neighboring gene, CYLD. We find associations with the IRF8 region and the region containing CDH1 and CDH3, as well as substantial phenotypic and genetic heterogeneity for CD itself. The genes are known to be involved in inflammation and immune dysregulation. These findings provide insight into the genetics of CD and suggest promising directions for understanding disease heterogeneity. The application of this method thus paves the way for understanding complex inheritance in general, leading to the dissection of different pathways and ultimately, personalized treatment.

PLAUR polymorphisms are associated with asthma, PLAUR levels, and lung function decline.

BACKGROUND: Several studies have suggested that chromosome 19q13.1-3 contains asthma susceptibility genes. OBJECTIVE: Linkage and association analyses using 587 United Kingdom and Dutch asthma families (n = 2819 subjects) were used to investigate this region. METHODS: A 3-phase procedure was used: (1) linkage and association analyses using 15 microsatellite markers spanning 14.4 mega base pairs (Mbps) on 19q13, (2) fine mapping of the refined region using 26 haplotype tagging single nucleotide polymorphisms (SNPs), and (3) dense gene analyses using 18 SNPs evaluated for association with asthma, bronchial hyperresponsiveness (BHR), FEV1, plasma urokinase plasminogen activator receptor (PLAUR), and rate of annual FEV1 decline in subjects with asthma. RESULTS: The microsatellite analyses provided tentative support for an asthma/lung function susceptibility locus (48.9-49.1Mbps), and fine mapping localized modest association to the PLAUR gene. PLAUR SNPs in the 5' region, intron 3, and 3' region are associated with asthma and BHR susceptibility and predict FEV1 and plasma PLAUR levels. SNPs in the 5' region showed association for asthma (2 populations), FEV1 (2 populations), and BHR (2 populations) phenotypes. SNPs in intron 3 showed association with asthma (2 populations) and BHR (3 populations). Importantly, the same 5' region and intron 3 SNPs were associated with plasma PLAUR levels. The same 5' region and 3' region SNPs were found to be determinants of FEV1 decline in subjects with asthma. CONCLUSION: This study represents the first report to identify PLAUR as a potential asthma susceptibility gene and determine PLAUR regions underlying this association, including a role in influencing plasma PLAUR levels. Finally, the association of PLAUR with lung function decline supports a role for PLAUR in airway remodeling in asthma.

Identification and replication of three novel myopia common susceptibility gene loci on chromosome 3q26 using linkage and linkage disequilibrium mapping.

Refractive error is a highly heritable quantitative trait responsible for considerable morbidity. Following an initial genome-wide linkage study using microsatellite markers, we confirmed evidence for linkage to chromosome 3q26 and then conducted fine-scale association mapping using high-resolution linkage disequilibrium unit (LDU) maps. We used a preliminary discovery marker set across the 30-Mb region with an average SNP density of 1 SNP/15 kb (Map 1). Map 1 was divided into 51 LDU windows and additional SNPs were genotyped for six regions (Map 2) that showed preliminary evidence of multi-marker association using composite likelihood. A total of 575 cases and controls selected from the tails of the trait distribution were genotyped for the discovery sample. Malecot model estimates indicate three loci with putative common functional variants centred on MFN1 (180,566 kb; 95% confidence interval 180,505-180, 655 kb), approximately 156 kb upstream from alternate-splicing SOX2OT (182,595 kb; 95% CI 182,533-182,688 kb) and PSARL (184,386 kb; 95% CI 184,356-184,411 kb), with the loci showing modest to strong evidence of association for the Map 2 discovery samples (p<10(-7), p<10(-10), and p = 0.01, respectively). Using an unselected independent sample of 1,430 individuals, results replicated for the MFN1 (p = 0.006), SOX2OT (p = 0.0002), and PSARL (p = 0.0005) gene regions. MFN1 and PSARL both interact with OPA1 to regulate mitochondrial fusion and the inhibition of mitochondrial-led apoptosis, respectively. That two mitochondrial regulatory processes in the retina are implicated in the aetiology of myopia is surprising and is likely to provide novel insight into the molecular genetic basis of common myopia.

Linkage disequilibrium maps and disease-association mapping.

Over the last few years, association mapping of disease genes has developed into one of the most dynamic research areas of human genetics. It focuses on identifying functional polymorphisms that predispose to complex diseases. Population-based approaches are concerned with exploiting linkage disequilibrium (LD) between single-nucleotide polymorphism (SNPs) and disease-predisposing loci. The utility of SNPs in association mapping is now well established and the interest in this field has been escalated by the discovery of millions of SNPs across the genome. This chapter reviews an association-mapping method that utilizes metric LD maps in LD units and employs a composite likelihood approach to combine information from all single SNP tests. It applies a model that incorporates a parameter for the location of the causal polymorphism. A proof-of-principle application of this method to a small region is given and its potential properties to large-scale datasets are discussed.

Effects of single SNPs, haplotypes, and whole-genome LD maps on accuracy of association mapping.

We describe an association mapping approach that utilizes linkage disequilibrium (LD) maps in LD units (LDU). This method uses composite likelihood to combine information from all single marker tests, and applies a model with a parameter for the location of the causal polymorphism. Previous analyses of the poor drug metabolizer phenotype provided evidence of the substantial utility of LDU maps for disease gene association mapping. Using LDU locations for the 27 single nucleotide polymorphisms (SNPs) flanking the CYP2D6 gene on chromosome 22, the most common functional polymorphism within the gene was located at 15 kb from its true location. Here, we examine the performance of this mapping approach by exploiting the high-density LDU map constructed from the HapMap data. Expressing the locations of the 27 SNPs in LDU from the HapMap LDU map, analysis yielded an estimated location that is only 0.3 kb away from the CYP2D6 gene. This supports the use of the high marker density HapMap-derived LDU map for association mapping even though it is derived from a much smaller number of individuals compared to the CYP2D6 sample. We also examine the performance of 2-SNP haplotypes. Using the same modelling procedures and composite likelihood as for single SNPs, the haplotype data provided much poorer localization compared to single SNP analysis. Haplotypes generate more autocorrelation through multiple inclusions of the same SNPs, which could inflate significance in association studies. The results of the present study demonstrate the great potential of the genome HapMap LDU maps for high-resolution mapping of complex phenotypes.

Genome scanning by composite likelihood.

Ambitious programs have recently been advocated or launched to create genomewide databases for meta-analysis of association between DNA markers and phenotypes of medical and/or social concern. A necessary but not sufficient condition for success in association mapping is that the data give accurate estimates of both genomic location and its standard error, which are provided for multifactorial phenotypes by composite likelihood. That class includes the Malecot model, which we here apply with an illustrative example. This preliminary analysis leads to five inferences: permutation of cases and controls provides a test of association free of autocorrelation; two hypotheses give similar estimates, but one is consistently more accurate; estimation of the false-discovery rate is extended to causal genes in a small proportion of regions; the minimal data for successful meta-analysis are inferred; and power is robust for all genomic factors except minor-allele frequency. An extension to meta-analysis is proposed. Other approaches to genome scanning and meta-analysis should, if possible, be similarly extended so that their operating characteristics can be compared.

Polymorphisms in a disintegrin and metalloprotease 33 (ADAM33) predict impaired early-life lung function.

RATIONALE: Asthma commonly originates in early life in association with impaired lung function, which tracks to adulthood. OBJECTIVES: Within the context of a prospective birth cohort study, we investigated the association between single nucleotide polymorphisms (SNPs) in a disintegrin and metalloprotease 33 (ADAM33) gene and early-life lung function. METHODS: Children were genotyped for 17 SNPs in ADAM33. Lung function at age 3 (n = 285) and 5 years (n = 470) was assessed using plethysmographic measurement of specific airway resistance (sRaw). At age 5, we also measured FEV(1). SNPs were analyzed individually using logistic regression, followed by linkage disequilibrium mapping to identify the causal locus. MAIN RESULTS: Carriers of the rare allele of F+1 SNP had reduced lung function at age 3 years (p = 0.003). When the recessive model was considered, four SNPs (F+1, S1, ST+5, V4) showed association with sRaw at age 5 years (p < 0.04). Using linkage disequilibrium mapping, we found evidence of a significant causal location between BC+1 and F1 SNPs, at the 5' end of the gene. Four SNPs were associated with lower FEV(1) (F+1, M+1, T1, and T2; p < or = 0.04). The risk of transient early wheezing more than doubled among children homozygous for the A allele of F+1 (odds ratio, 2.39; 95% confidence intervals, 1.18-4.86; p = 0.02), but there was no association between any SNP and allergic sensitization or physician-diagnosed asthma. CONCLUSIONS: Polymorphisms in ADAM33 predict impaired early-life lung function. The functionally relevant polymorphism is likely to be at the 5' end of the gene.

Positional cloning by linkage disequilibrium.

Recently, metric linkage disequilibrium (LD) maps that assign an LD unit (LDU) location for each marker have been developed (Maniatis et al. 2002). Here we present a multiple pairwise method for positional cloning by LD within a composite likelihood framework and investigate the operating characteristics of maps in physical units (kb) and LDU for two bodies of data (Daly et al. 2001; Jeffreys et al. 2001) on which current ideas of blocks are based. False-negative indications of a disease locus (type II error) were examined by selecting one single-nucleotide polymorphism (SNP) at a time as causal and taking its allelic count (0, 1, or 2, for the three genotypes) as a pseudophenotype, Y. By use of regression and correlation, association between every pseudophenotype and the allelic count of each SNP locus (X) was based on an adaptation of the Malecot model, which includes a parameter for location of the putative gene. By expressing locations in kb or LDU, greater power for localization was observed when the LDU map was fitted. The efficiency of the kb map, relative to the LDU map, to describe LD varied from a maximum of 0.87 to a minimum of 0.36, with a mean of 0.62. False-positive indications of a disease locus (type I error) were examined by simulating an unlinked causal SNP and the allele count was used as a pseudophenotype. The type I error was in good agreement with Wald's likelihood theorem for both metrics and all models that were tested. Unlike tests that select only the most significant marker, haplotype, or haploset, these methods are robust to large numbers of markers in a candidate region. Contrary to predictions from tagging SNPs that retain haplotype diversity, the sample with smaller size but greater SNP density gave less error. The locations of causal SNPs were estimated with the same precision in blocks and steps, suggesting that block definition may be less useful than anticipated for mapping a causal SNP. These results provide a guide to efficient positional cloning by SNPs and a benchmark against which the power of positional cloning by haplotype-based alternatives may be measured.